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Abstract. 

Recent investigations show that conservation laws limit the accuracy of gate 
operations in quantum computing. The inevitable error under the angular momentum 
conservation law has been evaluated so far for the CNOT, Hadamard, and NOT gates 
for spin 1/2 qubits, while the SWAP gate has no constraint. Here, we extend the 
above results to general single-qubit gates. We obtain an upper bound of the gate 
fidelity of arbitrary single-qubit gates implemented under arbitrary conservation laws, 
determined by the geometry of the conservation law and the gate operation on the 
Bloch sphere as well as the size of the ancilla. 
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1. Introduction 

One of the most demanding factors in realizing scalable quantum computers is the 
accuracy requirement for implementing the elementary quantum gates set by the 
threshold theorem for successful concatenated error- correct ion [1]. Two types of errors 
under current consideration are the environment-induced decoherence, caused by the 
interaction with the environment, and the controller-induced decoherence, caused by the 
interaction with the controller of the gate operation (the latter considered separately 
from, and in addition to, any errors arising from classical control imperfections). 
The environment- induced decoherence may be overcome, in principle, by developing 
qubits with long decoherence time. In the early treatments [1], the controller-induced 
decoherence was not distinguished from the environment-induced one. However, they 
have different constraints, since insensitivity and controllability are two contradictory 
demands. 

In the conventional description, the controller is described as a classical system and 
then causes no decoherence. However, the quantum nature of electromagnetic control 
fields has been studied by one of the present authors and others OIHlllllSlEllTllHllH], and 
shown to be a potentially substantial source of decoherence; see also Refs. [IDl [HI [12] 
for a debate on the validity of the model under consideration. 

In contrast to the above model-dependent approach, one of the present authors 
[T3| [HI [15], [16] independently found physical constraints on gate operations generally 
imposed by conservation laws, by quantitatively generalizing the so-called Wigner- 
Araki-Yanase theorem [T71 [T8l [T9] . The inevitable error probability under the angular 
momentum conservation law has been shown to be inversely proportional to the variance 
of the controller's conserved quantity for the CNOT gate [13] and the Hadamard gate 
|15] . while the SWAP gate obeys no constraint [13]. For the NOT gate a similar lower 
bound has been obtained more recently by a different method [20]. Subsequently, the 
above two approaches have been compared and merged [211 [22] . It was shown there 
that the limit derived from the angular momentum conservation law is equivalent to the 
one yielded by the phase fluctuations of the control field when a qubit-field interaction 
Hamiltonian of the Jaynes-Cummings type is assumed. 

Here, we extend the above results to arbitrary single-qubit gates with arbitrary 
conservation laws from a geometrical point of view. We show that a lower bound of the 
gate infidelity, one minus the squared gate fidelity, is given by 

sin^ f sin^ ^(1 - sin^ | sin^ 

4(l + a(L/c)2) ' 

where \1/ is the relative angle between the axes of rotations on the Bloch sphere generated 
by the conserved quantity and the gate operation; 6 is the rotation angle of the gate in 
the Bloch sphere, or the difference of the arguments of two eigenvalues of the gate, as 
defined later (in particular, 6 = ir for self-adjoint gates); L is the conserved quantity in 
the controller; and c is the maximum standard deviation of the conserved quantity in 
the qubit, while the standard deviation a{L/c) of L/c is generally understood as the size 
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of the controller. If the control interaction has full rotation symmetry, we eventually 
conclude that the gate infidelity of any implementation of an arbitrary gate on a spin 
1/2 qubit controlled by a spin N/2 ancilla system is bounded by (sin^6')(4 + 4A^^)~^ 
for < 6 < 7t/2 and (4 + AN^)~^ for 7r/2 < 9 < tt under the angular momentum 
conservation law. On the other hand, if the control interaction does not have the full 
rotational symmetry, the bound ([1]) holds. Because this bound vanishes in the important 
case 6 = IT and = 7r/2 (corresponding to a NOT, or bit-flip, gate, when the conserved 
quantity is proportional to the z component of the angular momentum), we have also 
developed here an approach that yields an alternative bound that is more suitable than 
([T]) under some circumstances; this alternative bound scales as (1 + ^/l + o"(L/c)^)~^ 
with the size of the controller. 

Although in many practical cases what we have called the controller will, by its very 
nature, be too large for our limits to represent a significant constraint, we note that, 
as pointed out above, laser pulses interacting with atomic qubits may be an exception 
to this [S], especially in view of the recent result that minimum energy pulses cannot 
be, in general, safely shared or reused [22]. Other systems for which our results may 
be quite relevant are the so-called "programmable quantum processors" [231 [2^ , which 
consist of a set of data qubits undergoing joint, closed evolution with a set of "program" 
qubits. If the total system evolution obeys a conservation law, our results clearly imply 
that the set of program qubits needs to be sufficiently large, in order to carry out the 
desired operations to a sufficiently large degree of accuracy. 

Our paper is organized as follows. Section 2 gives the formulation of the problem 
and examples to which our formulation applies. Section 3 introduces the deviation 
operator and gives a lower bound of the infidelity in terms of the variance of the deviation 
operator. In Section 4, we consider the commutator of the deviation operator and the 
conserved quantity to obtain a lower bound on the variance of the deviation operator 
using Robertson's inequality. In Section 5, the lower bound is represented in terms of 
the relative angle mentioned above and the size of the controller. Section 6 considers the 
case of the full rotational symmetry, while Section 7 gives an alternative lower bound 
which is useful for the case where the above bound vanishes. 

2. Basic formulation 

Let S be a qubit with computational basis {|0), |1)}. An arbitrary single-qubit gate 
is an arbitrary unitary operator Us on the Hilbert space Tis of S. The class of 
arbitrary single-qubit gates includes the Pauli operators defined by Xs = |0) (1| + 11) (0|, 
Ys = —i |0) (1| +i |1) (0|, and Zs = |0) (0| — |1) as well as the Hadamard gate defined 
by Hs = (v^)~^ (|0) (0| + |0) (1| + |1) (0| - |1) (1|). We introduce a vector a defined by 
a = {XsjYs, Zs) for convenience in later analysis. The vector a provides the general 
description of an arbitrary Hermitian operator of S 

Q 

A = <f)Is + -u-a, (2) 
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where m is a unit vector defined as m = {u^;, Uy, Uz) with WuW^ = + = 1, and 

(j) and 6 are real numbers, and Is stands for the identity operator of S. Here we have 
introduced an inner product between an arbitrary (c-number) vector m and the Pauh 
operators as 

u - a = u^Xs + UyYs + UzZs- (3) 

As an arbitrary unitary operator can be written as e*"^, with A Hermitian, we can write 
an arbitrary single-qubit gate f/g as 

Us = e^^ = e^'^ ^cos ^/s + i sin ij^ ' ^ ■ (4) 

Since two eigenvalues of Us are exp[2(0 ± 9/1)], parameters and 6 are uniquely 
determined with < < 27r and < ^ < tt. The angle 9 corresponds to the rotation 
angle in Bloch sphere induced by the gate Us] see p. 175 of Ref. p. 

Suppose we want to implement the arbitrary single-qubit gate Us hj letting the 
system S interact for a finite time interval with an ancilla A, which is described as a 
quantum system with a Hilbert space TYa- The total system to be considered is the 
composite system S + A, and we assume that its initial state is prepared as a product 
state before the interaction. We want the desired gate Us to be implemented by the time 
evolution of the composite system. Thus, every possible implementation is characterized 
by a pair a = (pa? U) consisting of a density operator pA on TYa describing the initial 
state of A and a unitary operator U on Hs ® Ha describing the time evolution of S + A 
during the interaction [13]. An implementation (characterized by) a = {pa,U) defines 
a trace-preserving quantum operation Sa by 

SM = TrA[U{ps®PA)U^] (5) 

for any density operator ps of the system S, where TrA stands for the partial trace over 
T^A- For simplicity, in what follows we will assume that pa is a pure state and 
in this case, the implementation will be described as a = {\A) ,U) (the case in which pA 
is an arbitrary mixed state is considered in the Appendix). 

How successful the implementation a = {\A),U) has been is measured by the gate 
fidelity of Sa relative to Us (Ref. [1], p. 418) defined by 

F{£a,Us) = miF{m, (6) 

W 

where \ip) varies over all state vectors of S, and F{\tp)) is the fidelity of the two states 
Usl'ip) and £ai\ip){'4j\), given by 

F{m = miSammus\i^y^'- (7) 

The implementation is perfect, i.e., Sa{p) = UspU^ for any density operator ps, if 
and only if F{Sa, Us) = 1. However, it has been shown that for several gates there is a 
constraint on the implementation naturally imposed by conservation laws [I3l [151 120] . 
In the present paper, we assume that there are additively conserved quantities Ls and 
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La of the systems S and A, respectively, so that the unitary operator U should satisfy 
the conservation law 

[U,L]=0, (8) 

where L = Ls ® Ia + Is ® -^a, and /a is the identity operator of A. 

Since different pairs of observables can represent the same additive conservation 
law, it is convenient to introduce a standard representation of the conserved quantity. 
Since a scalar operator poses no constraint, we always assume that Lg has two distinct 
real eigenvalues a < b. Let c = [b — a)/2. Then, it is easy to see that c is the maximum 
standard deviation of Lg attained, for instance, by \^) = 2~^/^(|Ls = a) + \Ls = b)). 
It is also easy to see that the operator c~^(Ls — bis) + /s is a non-scalar, unitary and 
self-adjoint operator. Therefore this operator can be represented as 

c-\Ls-bIs) + Is = l-(f (9) 
with a real vector / = (/^, ly, 1^) satisfying = 11 + 1^ + 11 = 1. Equivalently we have 

Ls = {b-c)Is + cl-a. (10) 

Here b and c are the maximum eigenvalue and the maximum standard deviation of 
Ls, respectively. The addition of a scalar operator to Ls does not affect the condition 
Eq. (IHl), and it also does not change the standard deviation of Ls; thus, the pair (Lg, La) 
and the pair [el ■ a, Lx) represent the same additive conservation law. From the above, 
it is also true that the pair (/ ■ a, La/c) represents the same conservation law as the pair 
(Ls, La). We shall use such simplifications where it is useful. 

As above, the conserved quantity Lg is determined by the vector /, which can be 
parameterized as / = (sin a; cos x, sin a; sin x, cos cj) with < u; < tt and < x < Stt. As 
a reference, we introduce a unit vector e = (0, 0, 1), which is the vector / for u; = x = 0. 
Then, there exists a unitary transformation Rs{l) of the system S such that 

Ls{l) = Rs{l)^ Ls{e)Rs{l) (11) 

for any /. For instance, we can take i?s(0 = exp[—'ico'(fi-(T)/2], where n = (/xe)/||/xe|| = 
(sin X, — cos X, 0) so that we have -Rs(0 = cos(u;/2)/s — z sin(ix'/2)(sinxAs — cosx^s)- 
Now, we assume that La also depends on I and that there is a unitary transformation 
Ra{1) such that 

La{I) = RAil)^Lj,{e)RA{l) (12) 
for any /. Then, the conserved quantity of the system S + A is transformed as 

L{r) = R(l)^L{e)R{l), 

where RiJ) = Rs{V) ® Ra{1) and L{V) = Ls{l) ® /a + /s ® ^a(0. 

The above formulation typically includes the following examples. 

Angular momentum conservation law. In this case, we assume that the system S 
has the spin {h/2)a. Then, Ls = {h/2)l ■ a represents the spin component along the 
(cu, x) direction. The angular momentum conservation law in the (u;, x) direction is 
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represented as the case where La is the angular momentum of A in the (cu, x) direction. 
In this case, we have b = c = h/2 and / is arbitrary. 

Atom-field interaction. An additive conservation law holds for the well-known 
Jaynes-Cummings model [25], which describes the coupling of a two- level atom S with 
a single-mode A of the electromagnetic field with annihilation operator a. Allowing for 
a detuning A between the atom and the field, the Hamiltonian for the model may be 
written in a suitable interaction picture as 

H = hAI (g) a^a + ihg (|0)(1| ® a - |1)(0| ® a^) , (13) 

where g is an appropriate coupling constant. Then, = Zg and La = 2a^a constitute 
a pair of additively conserved quantities for U = e~**^/'^ for any real t [21]. Thus, 
the constraint applies when one wants to realize a single-qubit gate by the Jaynes- 
Cummings interaction with the parameters A, g, and t. In this case, we have b = c = 1 
and / = (0, 0, 1). For multimode fields, we refer the reader to Refs. [211 122] . 

3. Mean square deviation 

In order to obtain an upper bound of the gate fidelity F{Sa,Us), we shall study the 
gate infidelity defined as 1 — F{£a,Us)'^, and show that an additive conservation law 
generally poses a lower bound on the gate infidelity. It is expected that such a limitation 
can be derived by the uncertainty relation generally formulated by Robertson [26] on 
quantum fluctuations, measured by standard deviations or variances, of arbitrary pairs 
of non-commuting observables. In fact, in Refs. [131 [15] commutation relations satisfied 
by noise and disturbance operators with the conserved quantity have been considered 
in order to apply Robertson's inequality, and lower bounds have been obtained for the 
gate infidelities of the CNOT gate and the Hadamard gate. Here, we extend the above 
method to arbitrary single-qubit gates. 

For this purpose, we introduce the deviation operator D of the system S + A defined 

by 

D = U\Ls ® Ia)U - Us^LsUs ® /a, (14) 

and we shall show that the variance of the deviation operator is, up to a constant 
factor, a lower bound on the gate infidelity. By the obvious cancellation when Eq. (fTOll 
is substituted in Eq. f|T^ . we can assume Ls = c / -c? without any loss of generality. Then, 
Lg has eigenvalues ±c. Let |xo) and be unit eigenvectors of Lg with eigenvalues c 
and —a, respectively. We define an orthonormal basis {l^o), l^i)} as |^o) = \xo) and 
l^i) = f/g Ixi)- For the input states with i = 0,1, the fidelity is given by 

Fm)= {X^\£amm\X^f'. (15) 

Since Tr[£^Q,(|,^j) (,^j|)] = 1, we obtain 

{X^\SA\Qm\X^) =l-F(|0))^ (16) 
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where i ^ j for i,j = 0, 1. The left-hand side can be described by using the general 
description of the output state f/(|^i) ® \A)) of the system S + A as 

U{\^i) ® \A)) = \xo) ® \Al) + \xi) ® \A\), (17) 



where \A^j) are unnormalized states of A. It is clear that 

{X^\SA\Qm\X^) = \\\Ai)r. (18) 

We, therefore, have 

iii^i)ir = i-^(ie,))' (19) 

for i ^ j. 

This fidelity F is related to the mean square (-D^) as follows. The mean square 
in the state ® \A) is written by 

(D^) = II [{Ls ® I)U - U{uILsUs ® mi^) ® \A))f. (20) 

Here, as any input state \ip) of S can be described as \ip) = cosCl'Co) + e*'^ sin ^ |,^i), 
where < C < f and < 6 < 27r, and Lg = c|xo)(xo| — c|xi)(xi|) from Eq. ([T7|l we 
have 

(Z}2) = 4c2( II I A?) f cos^ C + II \Al) r sin^ 0- (21) 
We substitute Eq. (fT9|) into this equation and obtain 

{D') = Ac' [I - F{\i,)f cos' C - i^d^i))' sin^ C] 

<Ac'[l-F{£^,Us)]. (22) 

Since <j{DY < {D'), where o-^D)' stands for the variance of D in the state ® 1^4), 
we obtain 

l-F{£^,UsY>^snpa{Df, (23) 
where varies over all the state vectors of S. 

4. Uncertainty relation 

In what follows, we shall consider the uncertainty relation between the deviation 
operator D and the conserved quantity L to obtain an lower bound of the variance 
of D. 

By the conservation law (IHl), we have 

[D,L]= -[Us^LsUs,Ls]®Ia. (24) 

Since addition of a scalar operator does not affect the commutator, we assume 6 = c or 
Ls = cl ■ a without any loss of generality. By Robertson's inequality [SB], we have 



aiD)aiL)>^\{[D,L])l (25) 
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where (■ ■ ■) and a stands for the mean value and standard deviation in the state 
respectively. Since the state {tp) ® \A) is a product state, the variance of the conserved 
quantity is given by 

a{Lf = a{Lsy + a{Lj,Y. (26) 

Thus, we have 

As the variance of Ls is upper bounded by the maximum eigenvalue c, namely, 
c"(-^s)^ < c, we see that 

„(Df > (1,) m^'^^"sM)\^ (28) 



v4c2 ; 1 + a(LA/c)2 

for any {tp) and \A). We are interested in the worst case error sup|^^ a{Dy'. Taking the 
supremum over \ip) of the both sides, and noting the relation sup|^^ |(-^)| = ||-^|| foi' 
any operator on Tis, we obtain 

SUpO-(L)) > — IT ; 77 -r--7i (29) 

for any \A). 

5. Constraint with relative angle in the general case 

In what follows, we determine the operator norm || [f/s^-^sf^s, -^s] || in geometrical terms 
from the vector analysis in three-dimensional space. 

To calculate this operator norm, recall the relations L^/c = I ■ a and Us = 
e"^(cos |/s + i sin |u ■ a), and calculate 

(f - a)(M ■ a) = (f - m)/s + i(lxiL)- a, (30) 

where I ■ u = IxUx + lyUy + IzUz and I x u = {lyUz — IzUy, IzUx — IxUz, IxUy — lyUx)- We are 
going to use Eq. (!30|) repeatedly. A straightforward calculation leads to 

U^LsUs = cos^ -Ls — 2c sin - cos - {I x u) ■ a + sin^ - {u ■ (r)Ls{u ■ a). (31) 

Now take the commutator of this with Lg. The first term gives zero, and the last one 
may be written as 

e 

2 I' 



sin^ -(u ■ a)Ls{u- a),Ls 



{I ■u){uxl)- {u x{lxu)}xl -a. (32) 



As u and / are unit vectors, the vector analysis gives {ux [l x u)) x I = —[I ■ u){u x I) . 
Thus we have 

[{u-a)Ls{u-a),Ls\ = Ac^i{l- u){u x I) ■ a, (33) 
and accordingly, 

[f/^Lsf/s,Ls] 
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4:ic^ sin 



cos ^ {(^ ^ X / 1 + sin ^ • (^u X fj 



a. 



(34) 



In Eq. f l34l) . the vectors (n x /) x / of the first term and {u x /) of the second term are 
mutually orthogonal and have the same norm \\uy< Hence the square of the operator 
norm of the commutator can be written by 

mlLsUs, LsW = 16cSin2 ^\\u x jcos^ ^ + sin^ ^ {)'} ' (35) 

because any vector a the operator norm \\a ■ a\\ is given by ||a||. 

Here, the vectors u and / have the following geometric meaning with respect to 
the Bloch sphere determined by a. We note that, in general, an operator of the form 
exp{—i6{n ■ a)/2) causes a rotation of the qubit state in the Bloch sphere by an angle 
6 around an axis given by the unit vector n; by the representation ([2]) and (jlj), we see 
that the desired gate operation Us is, up to a phase factor, precisely such a rotation, by 
an angle 6, around an axis given by u. On the other hand, the vector / represents the 
direction of the conserved quantity (for instance, the direction of a conserved angular 
momentum component), and therefore is an axis of symmetry of the system, since a 
rotation around /, given by Riit) = exp(— ztLs/2c), commutes with U. Let be 
the relative angle between those two axes corresponding to Ru and Rl. Then, \l/ is 
characterized by 

/■n = cos\l' and ||n x /|| = sin (36) 
Thus, we have 

II plLsUs, Ls]f = 16c^ sin^ ^ sin^ ^ (^1 - sin^ ^ sin^ . (37) 

We combine the above relation with Eqs. (!23l) and fl29|) . and obtain 

sin^|sin^^(l -sin^fsin^^) ^ ^ 

1 + (j{La/c)^ 



We have thus successfully generalized the result in Ref. [15] for the Hadamard gate to 
all the single-qubit gates. To see this, when 9 = n (corresponding to a self-adjoint gate 
Us), we have 

.0 sin^2\[' , , 

For Ls = Zs or Lg = Xs the relative angle of the Hadamard gate is 7r/4, and the lower 
bound becomes 4(1 + (t(La/c)2) \ which was proved in the paper [15]. When ^ = 0, 
the operator Us is proportional to the identity, and, of course, no restrictions apply in 
that case. 

Eq. fl37j) has a simple geometrical representation. See the left figure of Fig. [H In 
the Bloch sphere representation. Us is the rotation about the axis vector u by an angle 
9. By the unitary transformation of the conserved quantity, i.e., U^LsUs, the vector / is 
being rotated to the vector /' defined by U^LsUs = cV ■ a. Since || [a - a, 6- (x] || = 2 ||a x h\\ 
for any two self-adjoint operators a- a and h - a, we have || [f/gLsf/s, -^s] || = 2c^||/' x /||. 
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Figure 1. The left figure describes the rotation of the vector / in the Bloch sphere 
representation. The vector I rotates about the axis vector u by an angle 9, which 
corresponds to the unitary transformation uILsUs- The rotated vector I' and the 
original vector / make a relative angle 2j with each other. The right figure is just the 
interior elevation. 



Let the relative angle between the rotated vector /' and the vector / be 27, that is, 
||/' X /|| = sin 27. We rewrite the square of the operator norm as 

\\[uILsUs,LsW = Ac^sm''2^. (40) 

This equation must be consistent with Eq. fl37|) . This is easily seen from the right figure 
of Fig. [TJ Suppose that the point Q and R are located at the tips of the vector /' 
and /, respectively. The length of segment QR is described by sin 7, and at the same 
time it can be described by 2 sin sin | because PQ=PR=sin \E'. Therefore we see that 
sin^ 27 = 4sin2 | sin^ ^(1 - sin^ | sin^ '^), yielding Eq. ([37]) from Eq. 

In terms of the relative angle 27, the lower bound of the gate infidelity is written 

by 

l-F(g.,[/sf > ;^"7^^2V (41) 
4(1 + a{LA/c)^} 

If 6 = IT, it is easily seen from the left figure of Fig. [T] that 27 becomes 2\E'. In this case, 
we obtain the lower bound (139!) for any self-adjoint gate. This is expected because any 
self-adjoint gate is the rotation about the axis H by an angle vr. 



6. Constraint under the angular momentum conservation law 

Now, we consider the case in which all components of the angular momentum are 
conserved, which means that U must have full rotational invariance. In this case, for 
the given Us the conservation law ([HD holds for any direction {uj,x)^ so that holds 
for any ^ with c = h/2. We will derive a lower bound of the gate infidelity maximized 
over all the relative angles, which gives the lower bound under rotational symmetry. 
To get a lower bound, which holds for any initial state \A), we note that for any given 
ancilla system A, the standard deviation a{2LA^/h) is upper bounded by the operator 
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norm ||2La/^||, if it is finite. We obtain 

, sin^ fsin^ ^(1 - sin^ |sin^*) , , 

l + \\2Lj,{l)/hr 

for any initial state \A). Then the denominator of this lower bound is invariant over all 
the vectors /, because the operator norm is invariant under the unitary transformation 
defined by Eq. 

The numerator can be maximized by the proper choice of the vector / specifying 
the relative angle \Ef. Since the quantity sin^ | sin^ \E'(1 — sin^ | sin^ \E') is upper bounded 
by (sin^ e)/4 for < ^ < 7r/2, and by 1/4 for Ti/2<e < vr, we obtain 

l-F{£^,Usf> ^^^^X (43) 

4(1 + ||2La(/)//iP) 

for < ^ < 7r/2, and 

1 - F(^„ Usf > (44) 

4(1 + ||2La(/)//iP) 

for 7r/2 < < TT. The lower bound decreases as the operator norm of the conserved 
quantity increases. Thus, if A is a spin N/2 system we have ||2La/^|| = N and hence 

^ - nS.,U.f > (45) 
for < ^ < 7r/2, and 

l-F(£..,t/s)^>j(j^ (46) 

for 7r/2 < ^ < TT. In particular, if the ancilla consists of qubits, that is, the total 
spin number is N/2, we can conclude that any single-qubit gate on a spin 1/2 qubit 
system cannot be implemented within the gate infidelity less than (4 + 4A^^)~^ (up to 
the constant sin^ 6' for 0<6'<7r/2)bya rotationally invariant interaction with an N 
qubit ancilla system, or with any ancilla with spin N/2. 



7. An alternative lower bound 

The derivation in Section [5] (in particular, Eq. (l4Tll ) suggests that there is no constraint 
on the realization of gates for which \l/ = (that is, -u = /) or = 7r/2 (that is, 
u ■ Z = 0) and 6 = tt. While the first conclusion is correct, the second one is not true in 
general, since, for instance, a lower bound for the gate trace distance D{£a, Us) has been 
obtained recently [20] for the bit flip or quantum NOT gate. Us = Xs, under b = a = 1 
and / = (0, 0, 1). However, since the gate infidelity is always dominated by the gate trace 
distance, i.e., D{Sa, Us) > 1 — F{£a, Us)'^, we cannot immediately derive constraints for 
the gate fidelity from the result in [20]. Here, we shall instead show directly that such 
constraints do exist. 

To motivate the calculations that follow, consider again the case just mentioned, 
with Us = Xs (i.e., u = (1,0,0), 6 = ir, and (p = —tt/2), and a conservation law given 
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by Eq. fllOp with b = c = 1 and / = (0,0, 1) (i.e., Ls = Z^)- The difficulty is that in 
this case the ideal gate operator Us transforms Ls into something that commutes with 
itself: 

uILsUs = XsZsXs = -Zs = -Ls. (47) 

As a result, no constraint on D follows from Eq. ( l24|) . since the commutator on the 
right-hand side vanishes. However, the situation is different when we look at the action 
of Us on other system operators, in particular Xs and Yg. We have 

UgXsUs = Xs, 

UlYsUs = -Ys, (48) 
so if we define, as in Eq. dHj), the corresponding deviation operators 
D^ = U\Xs®Ia)U-Xs(»Ia, 

Dy = U^Ys ®Ia)U + Ys® /a, (49) 
we find (using the conservation law, [U, L] = [U'^ , L] = 0) 
[Dy, L] = U^Ys ® I A, L]U + [Ys, Ls] ® I a 
= 2iU^{Xs ®Ia)U + 2iXs ® I A 

= 2iD^ + 4iXs ® /a, (50) 

which shows that Dy and D^ cannot simultaneously be zero. Now, it is important 
to realize that the derivation of Eq. fl23|) does not depend on D being the deviation 
operator for the system's conserved quantity: an identical lower bound on the fidelity 
could be obtained from the deviation operator of any system Hermitian operator with 
eigenvalues ±c. Accordingly, in what follows we show how a generalization of the result 
( l50|l . to arbitrary Us and Ls, can be used in this way to derive a general lower bound 
for the gate infidelity that does not vanish when = 7r/2 and 9 = n, and thus in a 
sense complements the one obtained in Section 5, Eq. ( 14T|) . 

To simplify the calculations that follow, we introduce three operators /s,i defined 

as 

^s,3 = I ■ a, 
sm W 

, ... 1 r/r ^\ ri ^ 1 ^ ^ cos\E'- ^ 

/s,i = - ils,2ls,3 = -. — rl(f X n) X /} • cr = - — -u ■ a : — -/ ■ a. (51) 

sm W sm W sm W 

These operators obeys the same commutation relation as the Pauli operators: 

[ls,ijs,j] = 2i''^^eijkls,k, (52) 

k 

where eijk is the Levi-Civita symbol; in fact, they are essentially the Pauli operators, 
in a rotated reference frame. We define a new vector a' = {ls,iJs,2js,3)- Because any 
self-adjoint, unitary operator can be written in terms of a' as a • a', where a is a unit 
vector, we can write 
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with a unit vector u' = {ui,U2,u-i). From the definition (ICT) it is clear that 

U2 = ^Tr[/s,2(n' • a')] = ^-l^Tr[{(rx u) ■ a}{u ■ a)] = 0. (54) 
Therefore, Us may be written as 

Us = e*^ I cos ^Js + I sin ^-{u' ■ a') | (55) 

with U2 = 0. Straightforward calculations then yield 

Ulls,iUs = cos^ -ls,i - 2 sin - cos-{k x u') ■ a' + sin^ -(«' ■ a')ls,^{u' ■ a') (56) 

for i = 1,2, where we define li = (1,0,0) and I2 = (0,1,0). Here in the third term, 
{u' ■ ff')ls,i{u' ■ (?') can be written as 

(n' ■ a')ls,i{u ■ a) = (ni^ - U3^)ls,i + 2^1^3/8,3 (57) 

for i = 1, and 

{u-a')ls,2{u-a') = -{m' + us')ls,2 (58) 

for i = 2. We calculate the commutator [f/g/s,«f^S; ^8,3] by using the commutation 
relation (152|) and obtain 

( 6 9 9 9 ^ 

Vs,iUs, ls,3] = 2il- cos^ -ls,2 + 2u3 sin - cos -/s,i - {ui^ - us^) sin^ -/s,2 > , (59) 

f ^ 9 9 ^1 

[Usls,2Us, ls,3] = 2z < cos^ -/s,i + 2^3 sin - cos -/s,2 - {ui"^ + U's^) sin^ -/s,i > ■ 

(60) 

Let then Lg = (6 — c)/s + c/ ■ ? and define accordingly, by analogy with fj49l) . the 
deviation operators 

A = U^{ls,. ® /a)^/ - Ulls,,Us ® /a (61) 

for i = 1,2. Making use of the conservation law and the commutation relation (l52l) as 
in Eqs. ( I50ll above, we obtain the relation 

= -2zf/t(/s,2®/A)f/- [f/^/s,if/s,/s,3] ®/a. (62) 

As Ls can be written in terms of a' as Ls = {b — c)/s + c/' ■ a' with /' = (0, 0, 1), we have 
Ui = \\l' X n'll = sin^' and U3 = (/' ■ u') = cos\E'. These yield the following description 
of Eq. ([62D by substituting Eq. ([59]), 

^[Di, L] = -2i + 2Ulls,2Us ® /a + 2t7 • a' ® /a) , (63) 
where v = {vi,V2, fs) with 

Vi = 2 cos W sm - cos - , 
2 2' 

2 T • 2 ^ 2 ^ 

^2 = COS w sm - — COS -, 
9 9 

f 3 = — sin W sin - cos -. (64) 
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In a similar way, we have 

^[D2, L] = 2i (^Di + 2Ulls,iUs ® /a + 2w ■ a' ® /a) , 

where w = (^1,^2,^3) with 

2 T • 2 ^ 2 ^ 
wi = cos W sm - — cos -, 

e e 

W2 = — 2 cos W sin - cos -, 

T • 2^ 

= — sm W cos W sm -. 
^ 2 

Then, as in Section HI Robertson's uncertainty relation implies that: 
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(65) 



,2 ^ \{D2 + 2Ulls,2Us ^ Ia + 2v ■ a' ® Ia)\' 



(66) 



(67) 



1 + (t(La/c)2 

Now suppose we evaluate the numerator of the right-hand side in Eq. fl67|l in the 

eigenstate of Ulls,2Us corresponding to the eigenvalue +1. Let (■)' denote the 
expectation values in this particular state. We then have 

\2 + {D2y + 2{v ■ a')f 



Taking the square root of Eq. ( !68l) , we have 

As sup(X^) > a/ (X^) > holds for any observable X, we have 

/ j^j^^ ^ 2- v/sup((D2)^)-2|K7|| 

Vsup pi 2 > ,r / ' 

+ ct(La/c)2 

In a similar way, we have a lower bound of ^sup((-D2)^) as 

2- v/sup((Di)2) -2||«;|| 



VsUp((/^2)2) > 



v/l + a(LA/c)2 



(6^ 



(69) 



(70) 



(71) 



Adding Eqs. (|70D and (|7TD gives 



v/sup((Di)2) + v/sup((D2)2) > 



1 + v/1 + ^^(La/c)' 



:{4-2(||t;|| + ||.I;| 



(72) 

In the same way as in Section [3] (cf. Eq. fl22|) ). it can be shown that the gate 
infidelity is bounded from below by the mean square of Di, i.e.. 



l-F(^,,f/s)'>-^sup((A)'). 



4c2 



Therefore, we obtain 

l-F{8^,Usf> 



1 + v/i + ^(Wc)^ 



\w\ 



(73) 



(74) 
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Here I It; 1 1 and ||w|| can be written as 



9 



9 



9 



9 



cos^ — h cos^ \1/ sin^ - ( cos^ — h cos^ \1/ sin^ - 



9 



\w\ 



9 



9 



cos^ \1/ sin — h cos^ - ( cos^ — h cos^ \1/ sin - 
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(75) 
(76) 



This bound given by Eqs. (174l) - (176|) is fairly complicated. A slight simplification, 
resulting in a somewhat less tight bound, can be obtained by using 



\v\ + \w\<2 







\v\ H 


~ \w\ 



(77) 



This yields 

1 

> 

(l + ^l + a{Ljcf 

which can also be written as 




2 9 2 T 

1 — sin - sin W 
2 



2 - sin^ - sin^ ^ 
2 



(78) 



1 



1 + (r{Lj^/cy 



1 ^1 C0S7I 1 + cos^ 7 

v2 



(79) 



in terms of the angle 7 introduced in Figure 1. 

It does not seem possible to express the tighter bound given by Eqs. (17il) -( 176l) in 
terms of a single angle. However, in the case 9 = it, Eqs. f rM|l - fl76|) do yield a very simple 
bound: 



l-F{S^,Usf> 



1 



1 + ^1 + a{Ljc)'' 



1 



cos^|(l + I cos^l) 



.(80) 



The right-hand side of this equation is the lower bound for an arbitrary self-adjoint gate, 
i.e., f/g = Us (up to a phase factor); the Hadamard gate and the Pauli gates have this 
property. 

The general bound ( |7il) should be compared to the bound previously obtained by 
the first approach, that is, Eq. fl38|) . We plot these bounds both for 9 = n and 9 = it /2 
(See Figl2]). Generally speaking, neither bound is tighter than the other. For 9 = it/2 
and a{Lj^/c) = l,or 10, the previous bound ( l38l) is tighter than the present bound 
( !74|) over all the relative angle \E'. For 9 = it, the present bound does not vanish when 
= it/2; in fact, it reaches its peak, whereas the previous bound vanishes, which was 
the main motivation for this section. 



8. Concluding remarks 



We investigated the limitations of the gate infidelity of implementing arbitrary single- 
qubit gates under arbitrary additive conservation laws. We obtained the two different 
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Figure 2. Comparison between the lower bound in Eq. (fM]) (Solid line) and that in 
Eq. pS]) (dashed hne). We set 9 = n and (j{Lj^/c) = 1 in the upper left figure, and 
6 ~ TT and aiLp^jc) = 10 in the upper right figure, and d = 7r/2 and aiLxjc) = 1 in 
the lower left figure, and d — -k 12 and a^Lp^/c) — 10 in the lower right figure. 



lower bounds of the gate infidelity by using Robertson's uncertainty relation. Both 
bounds were described by the variance of the conserved quantity and the parameters 
{9,'^). The parameter \l/ is considered to be the relative angle between the axes 
of rotations on the Bloch sphere generated by the conserved quantity and the gate 
operation, which provides a geometrical understanding of the limitations in the three 
dimensional space. For the first bound (138|) . we show that it depends on the relative 
angle 27 between the axis of rotations generated by the conserved quantity before and 
after the gate operation, in the Heisenberg picture. Both bounds are important because 
in general, neither bound is tighter than the other. Both bounds become zero if and 
only if ^ = or \1/ = 0. This is expected because the conserved quantity L commutes 
with the desired gate Us only if ^ = or \1/ = 0. Generally speaking, the limitations 
decrease as the variance of the conserved quantity increases, and given the variance and 
the parameters 9 and the upper bound of both bounds provides the fundamental 
limitation for any implementation of the desired gate under the conservation law. 
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Appendix A. Ancilla in a mixed state 

In this section, for completeness, we show that our results are not changed if we assume 
an arbitrary state pa, which does not have to be a pure state, of the ancilla system as 
the initial state of the implementation Sa- In order to apply our previous treatment, we 
purify this state by introducing an auxiliary system B. The dimension of the auxiliary 
system B must be greater than or equal to that of the ancilla system so that the system 
B can provide the purification of any state of the ancilla system A. Suppose that \A') 
is the purified state for the state pa satisfying 

Trn[\A'){A'\] = pA, (A.l) 

where Tre stands for the partial trace over the auxiliary system B. \A') is a vector of 
the Hilbert space of the composite system A + B. This state can be considered as the 
initial state of the ancilla system if the ancilla system is extended to the original ancilla 
system A plus the auxiliary system B. In this case, the total system that we have is 
the composite system S + A + B, and the initial state of the implementation is given 
by ps ® ^he physical system is extended from the original system S + A to 

the system S + A + B to include the auxiliary system B, the time evolution operator 
of the implementation must also be extended to U' = U ® /b, where /b is the identity 
operator of B. Therefore the implementation is given by the trace-preserving quantum 
operation defined as 

Sa'ips) = TrA+B [U\ps ® \A'){A'\)U'^] (A.2) 

for any density operator ps of the system S, where a' stands for the characterization of 
this quantum operation, a' = {U', \A')), and TrA+B stands for the partial trace over the 
Hilbert space of A + B. The evolution operator U' must obey the conservation law 

[[/', L'] = 0, (A.3) 

where L' = Ls ® /a ® /b + -^s ® a ® /b j which is the natural extension of the conserved 
quantity according to the above extension of the physical system. This conservation law 
is equivalent to the previous one (IHl). 

It is not difficult to verify that this extension does not change the calculations 
leading to the general bounds, that is, we can apply the same arguments as from Sec. 3 
to Sec. 7 even in this case, and obtain 

sin^|sin^^(l -sin^fsin^^) , 
1 - F(g^,, Usf > ^ , , \, , J (A.4) 



and 



a(LA/c)' 

l-F{S^,,Usf> 



(l + + ^(La/c) 



^(IK'll + ll^l 



(A.5) 



The only difficulty here is that the implementation £a' appears different from the original 
one Sa- However, this problem is removed because the implementation Sa' is equivalent 
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to the quantum operation characterized by the time evolution operator U and arbitrary 
mixed states pa of the ancilla system, that is, 

Sa'ips) = TrA [U{ps ® Pa)U] (A.6) 
for any density operator ps- The proof is as follows. Since U' = U ® /b, we have 

SAps) = TrA+B [U ® /b(ps ® \^){^\)U^ ® /b] • (A.7) 
As |yl')(y4'| is the purification of the pa, we see that 

SAps) = TrA [f/TrB [ps ® U^] 
= TrA[f/(ps®TrB[|A')(A'|])f/t] 

= TrA[f/(ps®PA)f/], (A.8) 

which completes the proof. This proof has been already shown in Ref . [20] . 

We can then substitute Eq. (1A.6I) into Eq. ( 1A.4I) and ( 1A.5I) and obtain the general 



lower bounds, which are correct even if arbitrary initial states pa of the ancilla system 
are assumed: 

sin^fsin^^fl - sin^f sin^^) , 

1 - F{£^, Usf > ' fr / ^2 A.9 

1 + (t(La/c)^ 



and 



l-F{S^,Usf> 



1 



1-^(11^1 + 11^11)1', (A.IO) 



(i + v/i + ^(Wc)' 

where the quantum operation Sa is defined as 

£a{ps) = TrA [U{ps ® Pa)U] (A. 11) 

for any state ps- 
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